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Abstract This paper studies the positive solutions of a class of delay differential equations 
with two delays. These equations originate from the modeling of hematopoietic cell popu- 
lations. We give a sufficient condition on the initial function for t <0 such that the solution 
is positive for all time f > 0. The condition is "optimal". We also discuss the long time be- 
havior of these positive solutions through a dynamical system on the space of continuous 
functions. We give a characteristic description of the CO limit set of this dynamical system, 
which can provide informations about the long time behavior of positive solutions of the 
delay differential equation. 
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1 Introduction 

Delay differential equations are extensively used in modeling biological control systems, 
where the retardation usually originates from a maturation processes or finite signaling ve- 
locities t2ll8l l9lfTTIl . In this paper, we will consider the delay differential equation of form 



where T > Tq > 0, / 6 C(R+,]R+) and (j){t) > 0. This type of equation has been used to 
describe the dynamics of circulating blood cells 1 10 |. For obvious biological reasons, we are 
only interested in positive solutions of the equation. In this paper, we consider conditions 
on the initial function (j) (f ) such that the equation ([T} has positive solutions for all t > 0, and 
the long time behavior of these positive solutions. 
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- = -Yu + f{u{t - To)) - f{u{t - T))e-^(^-^o), t > 0, 



-T<f <0. 



(1) 



2 



Changjing Zhuge et al. 



The present work was motivated by investigating the dynamics of a mathematical model 
of hematopoiesis (4][5]|6|. This model consists of a set of four nonlinear delay differen- 
tial equations, describing the dynamics of proliferation and differentiation of hematopoi- 
etic stem cells, production of the three major types of circulating cells, leukocytes (white 
blood cells), erythrocytes (red blood cells) and platelets, and the feedback regulation of 
proliferation and differentiation ElU. These delay differential equations are obtained from 
age-structured models of each cell line using the method of characteristic line. For this set 
of equations, only positive solutions are biologically possible. In particular, the long time 
behavior of this system can be described by the co-limit set of these solutions, which is 
important for understanding the possible states of a system under given initial conditions. 

Bifurcations and bistability of the model of hematopoietic regulation have been studied 
numerically in ll3] |6llI0l . However, in numerical simulations we found that, in some cases, 
solutions of the delay differential equation system with a positive initial function can become 
negative at some time f > 0. This indicates that the delay differential equation model is not 
equivalent to the original age-structured model, which always yield positive solutions. The 
negative populations usually occur for erythrocytes and platelets, but not for stem cells and 
leukocytes. This is because the circulating erythrocytes and platelets are actively destroyed 
at a fixed time from their entering the circulating component |T'121. In out recent study 1 10 |, 
we proposed a set of initial functions such that the solutions are always positive. The current 
study provides a theoretical foundation for the findings in 1 10 ]. 

Equation (UJ can be obtained from a model of the population dynamics of erythrocytes 
or platelets |10|. In this model u{t) measures the population of circulating cells. The cell 
is produced from the differentiation of stem cells and amplification of precursor cells at a 
rate f{u(t)}. After differentiation, the cells undergo a stage of maturation of duration Tq, 
and then enter into the circulation. The circulating cells are lost at a rate y, and are actively 
destroyed at a fixed time (t — Tq) from their time of entry the circulating compartment. 

This paper is organized as follows. In section |2l we prove a sufficient condition for 
the existence of positive solutions of (UJ. At the end of this section, we will show that the 
condition is "optimal" by an example. In section[3] we will introduce an iterated map on the 
space of continuous functions and study the co-limit set of positive solutions of (O through 
this map. The paper concludes with an example in section|4] 

2 Sufficient condition for positive solutions 

In this section, we develop a sufficient condition for the initial function (j){t) such that the 
solution of equation (Ujl is positive for all time f > 0. 

Throughout this paper, we define an operator H : C([— T, +o°),R+) x ]R+ — ^ R such that 

H(u,t) = u(t)- r f{u{t-a))e-^'-''-''°Ua, (2) 

and let//o : C([-T, +°°),K+) ^ M be defined as 

Hq{u) = H{u,0) = m(0) - r f{u{-a))e-T^'^'-'">Ua. (3) 

Theorem 1 Consider the equation (|T]l. If the initial condition (j){t) satisfies 

Ho{^)>0, (4) 
then (O has a unique solution u(t) in t £ [— T,oo), and the solution is positive for all t > 0. 
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Proof First, the existence and uniqueness of the solution is straightforward using the method 
of steps 

Next, we will show that the solution u{t) of equation (UJ can be given iteratively through 

u{t) = e-^HQ{i<)+ r f{u{t-a))e-^^''-^°Ua. (5) 

J To 

From (|T]|, it is easy to obtain the following when t >0, 



4t) = e-^ (u{Q) + ^' [/(K(a - To)) - f{u{a - t))^-'''"-"'')] e^"da 



(6) 



Thus, we have 



g7(fl-T+To) 



M(f) = g-r' /^„(0) + f{u{a-x))e^^''-^^^^da- f f{u{a- T:))e^'-''-^+^'>Ua] 

\ Jx-xq Jo J 

= e-^'(u{0)+r f{u{a-T))e^^"-'+'°^da+ [' f{u{a- t)) 
= e-^ (u{0)- r f{u{-a))e^^-''+^o)aa+ f f{u{t-a))e'^^'-''+^''Ua 

J To 

Here, the initial condition u{t) = (j){t) has been applied in the last equality. 

Now, we have derived ([5J, from which the solution is always positive for all f > provide 
Ho{(j)) > 0. Thus the theorem is proved. □ 

Remark 1 Theorem [T] can be extended to the case of time-dependent function / = f{u,t). 
In this case, equation ([TJ becomes 

^ = -YU + f{u{t - To), t - t:o) - f{u{t - T:),t - T)e-^'^^-'''>\ t > 0, 

u{t) = (j){t), -T<f<0. 

Accordingly, the operator H is redefined as 

rT 

H{u,t) = u{t)- f{u{t-a),t- a)e'^'''''^Ua, 

J To 

andHo{u)=H{u,0). 

Similar to the proof of Theorem [T] we can show that the solution of Q is given by 

u{t) = e-^' Ho{u) + r f{u{t-a),t-a)e-^^''-^»^da. (8) 

Jtq 

Thus, if HQ{(j)) > 0, the solution of (|7]l is positive for all 1 > 0. 

Remark 2 The condition © is not necessary for a given function /. However, we argue that 
this condition is 'optimal' for general /, i.e. for any e > 0, there exist a pair of functions, / 
and (j), such that 

Ho{<j)) < -e, 

and the solution of ([T]l with initial condition (j){t) becomes at some r* > 0. A construction 
of / and (j) is given in Proposition[T]below. 
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Proposition 1 For any e > 0, let (^{t) = (po > {t e [— T,0]), and 

/a(«)=Ae-("-*«)'/^ (9) 

where 5 > and 

^("^"t^) . (10) 

Let ug (t) be the solution of 

^ = -r« + /6(«(r-To))-/5("(?-T))e-''(^-^«\ f>0, ^^j^ 
u{t) = (j){t), -T<r<0. 

Then when 5 is small enough, the solution satisfies ug{T) < 0. 

Proof When 5 — > 0, the function fg (m) converges to the following function 

Let uo{t) be the solution of jllb with fg replaced by /q. 

We will divide the proof into two steps. First we will prove that 

lim ug{t) = Mo(f) (13) 

for all t > 0. Next, we will show that uq{z) < 0. Therefore, we have ug{T) < when 5 is 
sufficiently small. 

To prove l ll3b , we first note that the proof of Theorem [1] is also valid when f{u) is 
replaced by /o(m). Thus, «()(/) also satisfies the iterative equation Therefore, for any 
? > 0, we have 

lim \ug{t) - UQ{t)\ 
= lim e-r- (^£[/5(<^(-«))-/o(<^(-fl))],-r(«-xo)j«^ 

+ r[fg{ug{t-a))-fo{uo{t-a))]e-'^<^''-''>Ua . 
< lim / \fg{^Q)- fQ{^Q)\da+\\m I \fg{ug{t -a)) - fQ{uQ{t -a))\da 
= f l^r[i\fg{ug{t-a})-fo{uQ{t-a})\da. 

In the last step, we note that \ fg{ug (t — a))— fQ{uo{t — a)) \ < 2A, and therefore the equality 
holds according to the Lebesgue dominated convergence theorem. Now, we will prove i ll3b 
inductively. 

When < f < To and Tq < a < T, we have ug [t — a) = mo(? — a) = (^o- Therefore 

lim\ug{t)-uo{t)\= [ lim \fg{ug{t -a}) - fo{uo{t -a})\da = 0. 
5-s-O .JtQ 5-s-O 
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Now, we assume that ( 113b is valid for < f < wTq. When wTq <t< {n+ 1)tq and Tq < 
a < T, we have (f — a) < iiTq, and therefore 



lim us{t — a) = UQ{t — a) , wTq </<(n+l)To, To<fl<T. 

5->0 



Thus, we have 

lim |Mg(f) — Mo(f)| = / \\m\f§{u§{t — a)) — fQ{u(i{t — a))\da = Q, 
5->0 Jtq 5-s-O 

and ( IT3l l holds for < f < (n + 1)to. Therefore, we have ( IT3] l for all ? > 0. 
To prove that mo(t) < 0, we will show that 

f f{uQ{t-a))e-^^''-^^Ua = Q. 

JtQ 

In fact, we will prove a stronger result. Let 

/„ = [«To,(n+l)To], n = 0,l,2,---[T/To] + l, 

where [t/tq] denotes the integer part of t/tq. We claim that in each interval /„, there is at 
most one point f„ £ /„ such that uo(t„) = (j)o. 
First, we note that dlOl l is equivalent to 

Ho{u) = (1)0- f f{^o)e-y^''-'^Ua = -e. (14) 

7 To 

From ([5]l, when < ? < Tq, 

JtQ 

J to 

= (^)o + e(l-e-^)>(?o. 

Thus, u{t) (j)o for any t 6 /q, and hence, 

/(Mo(r)) = 0,Vfe(0,To]. 

Now let n < [t/tq] + 1, and assume that for any integer k < n, there is at most one point 
ti; E 4 such that Mo(fji:) = <?o- For f 6 /„ and f < T, we have 

UQ{t) = -ee-^+ r f{uo{t-a))e-^^"-'»Ua 

J To 

We note that when t 6 7„ and Tq < a < T, < f — a < wTq. From our assumption, there are at 
most n points in f 6 [0,wTo] such that uo{t) = (j)Q. Therefore, the first integral 

[' f{uQ{t-a})e-^^"-">'>da = 0. 

J To 
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Thus, 



(0 = -ee-^' + f{uQ{t-a))e-^^"-^^Ua. 



Hence, 



^ = e-r'(ye -Ae^^o) < -e-^'Y(t>o < 0, 



and uo{t) is strictly monotonically decreasing for t E I„. Therefore, there is at most one point 
t„ e /„, f < T such that u{t„) = and the claim is proved. 

Now, since there are at most a finite number of values of ? 6 [0, t] such that uo{t} = (j)Q, 
the integral 



/V("o(T-fl))e-''("-"°' 

7 To 



da = 0. 



Therefore, 



uo{t) = -£e-^^+ r f{uo{t-a))e-^'^"-''«Ua=-ee-^^ <0, 

J to 

and the proposition is proved. 

Figure [T] shows the solution of jllb when 5 is sufficiently small. 




Fig. 1 An example solution of the equation il 11 . Parameters used are 7=0o = l>To=O.l,T=l,e = l,5 = 
0.005. 
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3 (B-limit set of positive solutions 

In this section, we always assume that / is bounded, 

0<f{u)<f, VMe]R+, (15) 

and discuss the co-limit set of positive solutions. 

It is easy to see that any solution u(t) of Q is associated with a sequence of functions 
{M„(f)}„gN in C([— T,0],K+) such that u„{t) = u{t + nT). From ([5]l, we define a map T : 
C([-T,0],K+) C([-T,0],]R+) such that m„ = Tun-i- Explicitly, T is given by 



{Tu){t) 



'^Hq{u)+ r f{u{t + T-a))e-^^''-'»'>da, fe[-T,To-T], 



J To 

+ [ /(M(f + T-fl))e-^("-^»'da, fe(To-T,0] 
Jt+i: 



(16) 



where 

Ho 



(u) = u{0) - r f{u{-a))e-T^''-^»'>da 

J to 



as defined in ((Sjl. Therefore, T defines a dynamical system in the function space C( [— T, 0] , 
For any mq 6 C([— T,0],]R+), the co-limit set of uq under the map T is defined as 



C0r(w()) = < V 6 C([— T,0],R ) There is a sequence — > -1-°° such that lim T"''uq = v 

The structure of the co-limit set c07-(mo) provides characteristic descriptions of the long time 
behavior of the solution of (|T} with a given initial function ii{t) = uo{t) when — T < r < 0. 

According to the previous discussions, we define two subsets of C([— T,0],K+) as fol- 
lows. 

^ = {MeC([-T,0],R+) |//o(m) >0}, (17) 
% = {ueC{[-T,0],R+)\Ho{u)=0}. (18) 

Further, we equip & and &o with the norm: 

||m|| = inax 
ie[-T,0] 

It is obvious that &o C i^. 

Theorem 2 Let T be the operator defined in ( 1161 , we have 

1. T{&) C ^. 

2. T is completely continuous in the C([— T,0],R"'") topology. 

3. For any u (z 3!, {r"M}„gp^ is uniformly bounded. 

4. For any m € the (0-limiting set (Or («) is nonempty. 

5. For any u e COt{u) C 

To prove Theorem|2l we need several lemmas. 
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Lemma 1 For any fixed t 6 M^, the operator H(u,t) is continuous with respect to u on the 
Banach space (C([— T,+o°),IR+), || • ||) with the norm 



sup \u{t)\. 

te\-x.+oo) 



In particularly, Hq{u) is continuous. 



Proof Let mq be any function in C([— T,+°°),]R+). Since uq(s) is continuous for s > — T. 
For any f > 0, minTj,<fl<T UQ{t — a) and max^^<^a<x "o(' ~ 'j) are well defined and finite. 

For any e > 0, since f{u) is a continuous function, it is uniformly continuous when u 
takes value from the interval 

£ £ 

min u<\{t — a) , max ui\(t — a)-\ — 

To<a<T 2 To<a<T 2 

Thus, there exists 8 : < 5 < £/2, such that for any u e C([-T, +°°),IR+) with ||m-mo|| < 8, 
and hence u(t — a)El for any a 6 [tq, t], have 



\f{uo{t-a))-f{u{t-a)) \ < 



2(t-To) 



Va e [To,t]. 



Thus, 



< 



H(u,t)-H{uo,t)\ 

u{t) - uo{t) - r [f{u{t -a))- f{uo{t - a))] e-^^''-'o^da 

J To 

||m — Mo||+ / \f{u{t — a))—f{u(){t — a))\da 

Jtn 

da = e, 



£ 



£ 



2 Ao2(t-To) 
and the lemma has been proved. 

Lemma 2 For any u{t) e 

Ho{T"u) = e-^"^Ho{u), Vn e N. 
Proof From the definition of T in i ll6t . we have for any u E 

{T"u){0)=e-^'Ho{T"-\)- r f{{T"u){-a))e-^^"-^'>'>da. 

JtQ 

Hence, 



(19) 



Ho 



{T"u) = {T"u){0)- r f{T"u{-a))e-^^''-''°' 

J To 



da 



Thus, il% is obtained by induction. □ 

Lemma 3 For any uo{t) € & and the sequence {T"uo}„^ff, let the function u{t) on [0, 
be defined by 

u{t) = {T"uo){t-nT), t e [(n- 1)t,wt),« £ N. 
Then u[t) is the solution of (O with initial function u(t) = mo(0 when — T < / < 0. 
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Proof First, we show that u(t) is continuous. To this end, since u{t) is continuous on each 
interval ((n — l)T,nT), we only need to show that u{t) is continuous at all points t = nz . For 
any v 6 from il6l . we have 

(rv)(-T) = v(0) - r f{v{-a))e-^^"-''>'>da+ f f(v{-a))e-^^"-^»'>da = v(0). 



Therefore, 

lim M(o = (r+^(o)(-T) = (rMo)(o)= lim m(0 = (r+iMo)(-T), 

and is continuous at f = «T. 

Next, we will prove that u{t) satisfies equation l|5ll: 



Jtn 



(20) 



To this end, we only show that ( 120b holds for (« — 1 ) T < f < (« — 1 ) T + Tq (Vw 6 N) , the case 
when (« — 1)t + To < r < nT is similar. 

When {n — 1)t < t < (n — 1)t + To, we note that for a e [to, t], (n — 2)t < f - a < 
(« — 1)t, and therefore 

{T"'\o){t-a{n-l)T) = u{t-a). 
Thus, applying Lemma|2l we have 

u{t) = (r(r-iMo))(f-«T) 

= e-''('-"^+^)//o(r""'"o)+ f f{{T"-\io){t-m + X-a))e-'^''''-^Ua 



= e-7('-(«-i)^)e-7(«-i)^//g(„g)+ /'V((7'""'Mo)(r-a-(«-l)T))e-''(''-"'''^/a 
= e-'"i/o(Mo)+ r f{u{t-a))e-^^"-^'>Ua, 

J to 

and l|20]l holds. 

Now, we will prove that u{t) is the solution of ([U with the initial condition Mo(f). It is 
obvious that u{t) = mo(?) when — T < f < 0. When r > 0, from i l20b , we have 



dt dt 

d f'-'^o 



-Y{t-to} 



= -ye-'^'Hoiuo) + — [ ° f{u{a))e-^'''-''-'»Ua 
dt Jt-x 

= -ye-^' Ho{uo)+f{u{t- To)) - f{u{t - T))e-^^'-'°^ - Y r'"fHa))e-^^'-^^^^^ 

Ji-t 

= -Y(^e-r'Ho{uo) + J^' f{u{t-a))e-r^"-'"''da^ +/(«(^ _ Tq)) -/(«(?- T))e-^( 
= -Yu{t) + f{u{t - To)) - f{u{t - T))e-^(^-^«' . 

Thus u{t) is a solution of l[TJ. 

Lemma 4 For any u (z {Tu){t) is continuously differentiable. 



10 



Changjing Zhuge et al. 



Proof For any m 6 i^, it is easy to see that {Tu){t) is continuous in f 6 [— T,0]. From the 
proof of Lemma [3] we have 



-YiTl,){t)+fiu{l + T-To))-f{u{t))eyi'-'o\ te[-T,To-T], 

(21) 

-Y{Tu){t) + f {{Tu){l + T- To)) -f{u{t})eri'-'o),ie {To -T,0]- 
Therefore, Tu is continuously differentiable. □ 
Now we are ready to prove Theorem (2] 

Proof (Proof of Theorem^ 

(1) For any m 6 i^, it is easy to have Tu e C([— T,0],IR+). In addition, from ( 116b , we 
have 

{Tu){0) =Ho{u)+ r f{{Tu){-a))e-^'^"-'»Ua 
J la 

> r f({Tu){-a))e-^^''-^»'>da>0. 

J to 

Thus, Tu e and hence T& C 

(2) To prove that T is completely continuous, we only need to prove that for any bounded 
sequence {m„}„>i C {Tu„}n>i has a convergent subsequence. To this end, from the 
Arzela-Ascoli theorem, we only need to show that the set of continuous functions {Tu„}n>i 
is uniformly bounded and equicontinuous. 

Assume that ||m„|| < M < o° for any n > I. From the fact that / is bounded, i.e. the 
inequality ( 115b , and the definition of T in ( 116b , we have 

\Tun{t)\ < |m„(0)|+2 r fe-^^"-^°Ua, \/t € [-T,0]. 
Jta 

Thus, 

||rM„|| < |1m„||+2/(t-To) <M + 2/(T-To). 
which implies that {Tu„]„>\ is uniformly bounded. 

The sequence {Tu^} is equicontinuous if the sequence | — {Tun) 1 is uniformly bounded 



in f 6 [— T,0], which is proved below. From ( 121b in Lemma|4l it is easy to have that 



d 

— TUnit) 



< r\\Tu„it)\\ + If < Y{M + 2/(T - To)) + 2/ 



is bounded. Thus, {Tu„} is uniformly bounded and equicontinuous, and therefore T is com- 
pletely continuous. 

(3) From ( 116b . for any m 6 i^, we have 

\{T"u){t)\ = \T{T"-'u){t)\<e-r^'+'^Ho{T"-'u) + {T-To)f, Vre[-T,0]. 

From LemmaEl Ho{T"-\) = e-^^"-^'i''Ho{u), thus, 

||(r"M)||<//o(M) + (T-To)/, (22) 

which implies that {T"u]„i^^ is uniformly bounded. 
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(4) Since the sequence {T''u}„(^n is bounded in C([— T,0],IR+) and T is completely 
continuous, there exists a subsequence rik — ^ +°° such that T"''u converges to a function in 
C([— T,0],K+). Hence for any u E 0}r{u) is nonempty. 

(5) For any v E COt{u), there exists a sequence nj^ — > +oo such that v = lim,t^+oo r"*M. 
From Lemma[T] Hq{u) is continuous with respect to u. Thus, we have 

Ho{v)= Mm Ho{T"''u)= Mm e-^"'<''Ho{u) =0, (23) 

and hence v 6 i^o- Therefore (Oriu) C £^o- D 

From Theorem |2] informations about the long time behavior of positive solutions of (UJ 
are contained in f^o- In particular, all solutions of ([TJ with initial function in will be of 
great interest. In rest of this section, we will focus our discussion on ^o- 

From Lemma |2l we immediately have the following corollary. 

Corollary 1 The restriction of T on &() is an operator from into &q. 

Proof From Lemma[2l for any u 6 !^q, we have Hq{Tu) = e^^'^HQ{u) = 0, and hence, Tu e 
%■ □ 

Corollary [T] implies that if m £ i.e. Hq{u) = H{u,0) = 0, then H{u,n'c) = for any 
n £ N. In the following Proposition |2l we prove further that if H{u,t*) = at one time point 
;* > 0, then H(u,t) = for all r 6 R+. 

Proposition 2 Assume that u(t) is a solution of lO, then 

H{u,tx)=e-^^'^-'^'^H{u,t2)- (24) 
In particular, if there is f* > such that 

H{u,t*) = Q, (25) 

then ( 1251 is satisfied for any t 6 . 
Proof From ([5J, we have 

H{u,t) = u(t) - r f{u{t-a))e-^'^"-''"^da = e-^'Ho{(t>) = e-^'Ho{u). 

Thus we have 

//(M,fi) =e-^<'i-'^'//(M,r2), Vri,f2 > 0. 
The second part is obvious when we take t2 = t* in l l24b . □ 



Remark 3 From i l24b . we have 

F{u,t) = e^'H{u,t)=Ho{u) (26) 
for all t > 0. Thus, F{u,t) is a 'first integral' of the delay differential equation ([T]l. 
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Remark 4 In i24i . if we take ti =t, and t2 = 0, then 

H{u,t)=e-^Ho{u) 

which implies 



dt 



{H{u,t)) = -YH{u,t). 



We note that Hq{u) > for any m e f^. Thus, H{it,t) is a 'Lyapunov functional' of the 
equation ([U on Furthermore, for any solution of (|T]| with initial function (j) E S>, the 
solution approaches &o on a time scale 1 /y. 

In the proof of Theorem[2l we have shown that is nonempty. The simplest case is for 
a constant function u{t) = uq. Then uq should satisfy the equation 



MO ■ 



du 
dt 



: rf{uo)e-^^"-'''Ua = f{uo){l-e-^'')/Y, 

J to 

i.e., M = Mo is a steady state of the equation 

+ _ To)) - /(M(f - T) )e-^(^-^») . (27) 

It is not trivial to find a non-constant function in &o, i.e. a function u{t) such that 

m(0)= r f{u{-a))e-^''''-''oUa. (28) 

Proposition [3] gives an explicit way to construct a function in ^q, which is important when 
we solve equation |T]l numerically and perform a bifurcation analysis. 

Let g{t) be a function in C([0, t],]R+). We define an operator Fg : C([0,t],M+) 
C([-T,0],K) as follows. 

When Ti = T — To > To, F„M is defined as 



fiu(t- To)) -g{-t)e-y'^, -To<r<0, 
{Fgu){t) = { f(^u{t-To))-g{-t)e-^^'+'\ -Ti<r<-To, 
g{-{t + Ti))-g{-t)e-y('+'\ -T</<-Ti. 
When Ti = T — To < To, FgM is defined as 



(29) 



{Fgu){t) = { 



f{u{t-To))-g{-t)e-^'^, 
g{-{t + ri))-g{-t)e-r^^, 



-Ti < r < 0, 
-To <r < -Ti, 



(30) 



?{-{t + Ti))-g{-t)e-y('+'\ -T<r<-Ti. 



It is easy to see that {Fgu){t) is a continuous function of t in [— T,0]. 
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Proposition 3 For any g e C([0, t],IR+), let u{t) ( — T <t <Q) be the solution of the differ- 
ential equation 

du 
dt 



^(-T) = f g{t)e-y'dt. 

J To 

Then u(t) € If further, g satisfies 



— = —yu + FgU, 



(31) 



g{0)=f(^J\{t)e-y'dt^ 



(32) 



then, u{t) is continuously differentiable, and 



lim = -Yu{0)+f{u{-To))-f{u{-T))e-^^'-'«\ (33) 

Proof We only prove this for the case Ti > Tq, the Ti < Tq case is similar. 

We just need to solve the equation ( 13 lb directly and verify ( 128b . Using variation of 
parameters, we have 



u{t) = e-^ ( r g{t)e-T^^dt+ f e^'{Fgu){s)d! 



Letting f = 0, we have 

i(0) = f g{t)e-'"dt+ r\v^g^_, + ^^)_g^_,),-7[s+r)y^ 



u(V} = I g(t}e '-at-t I ' e'-(g(-s + T:i}-g(-s}e ''-"'jas 

To 

+ / '\-^{f{u{s-To))-g{-s)e-''^'+'^)ds 



r\r'{f{u{s-To))-g{-s)e 
f er'{f{u{s-TQ))-g{-s)e-^'')ds 
f\{t)e-^'dt+ ( '\{-s + Tx)e'"dt- [ '\{-s)e-^'ds 
+ f f(u[s-Xo))e'"ds- f g{-s)e^^'-^^Us 
rg{t)e-^'dt+ f g(-r)e^<''-^i)jf- / '\{-t)e-^'dt 
f f(u{s-Xa))e'''dt- f g{-t)eT^'^'-''Ut. 



+ 



= r g{t)e'^^dt- r g{s)e-T^^dt+ r f{u{-t))e-^^'-''°Ut 

J Zq J Zq J Zq 

= r f{u{-s))e-^^'-'»'>dt. 

Thus, u E &o and the first part is proved. 

rr, , . . r du , , , du . 

To prove the contmuity oi °nly need to show that — is contmuous at f = —X\. 

From ( 132b . 



g(Q)=f(^j\{t)e-y'dt^ =/(«(-T)), 
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we have 

du J . , . 

Thus, — is continuous at f = — Ti from Olb and ( 129b . 

at 

To verify i33i . we note that 

hm u'(t) = lim ( -Yu{t) + f{u{t - Tq)) - g(-t)e^^''^) 

= -Yu{0)+f{u{-To))-g{0)e-r'^ 

= -ruiO)+fiui-To))-f{u{-T))e-^^'-'»l 

The theorem has been proved. □ 

Biologically, the function g{a) in Proposition[3]corresponds to the initial age distribution 
of cells differentiation from stem cells. From Propostion[3] if we choose a function g{a) that 
satisfies ( 132b . then the solution of ( 131b can serve as an initial function of the equation ([TJ, 
and the solution of (O therewith is positive for all f > 0, differentiable in f 6 (— T,+oo) and 
satisfies H{u,t) = for all t > 0. 



4 Examples 



In this section, we will show numerical results for two examples, with 

/(m) = ^— , (34) 

and 

2 

/(m) = — arctan(asin(2M:M)) + 1, (35) 

71 

respectively. The co-limit set is illustrated at figure |4] in which N{cOt) is a map from M to 
subsets of K defined as 

N{o>r)= U U {[°^if)dt\. 



In the simulation, we choose m e i^o randomly according to Theorem [3] and generate a 
sequence {T"u} with n large enough (we calculate to n = 600) till the sequence converges. 
We take the last 100 functions from {T"co} as the ©-limit set C0t{u}. For each v 6 COt{u), 
we calculate the L' norm ||v||^i = J^^v{t)dt and then give one point in Figure|4] 

Figure |4] can be thought of as a 'bifurcation diagram' of the delay differential equation 
(O. When N{cOt) has only one point, it means that the system has a globally stable steady 
state. When N{(j}t) contains a finite number of points, then either the system has multiple 
steady states, or there is a periodic solution, with period qT for some q a rational number. 
When N{cOt) contains an infinite number of points, then the system either has a periodic 
solution with period qz and q is an irrational number, or chaotic solutions. Chaotic solutions 
have been found in many simple delay differential equations f7l fTT]|I3l , including the well 
known Mackey-Glass equations, which correspond to the case with T = +°° in the present 
study. The approach in this study provide a way to study complex behaviors of such systems 
from the point of view of dynamical systems in a function space. 
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Fig. 2 Illustration of the a)-limit sets of positive solutions of (T) with nonlinear function / defined as (34) 
and i35\ . respectively. Black points are L' norms of stable steady states or oscillatory solutions, red points 
for unstable steady states. Parameters used are To=0.1,T=l,7 = 0.005. 
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